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We present a measurement of the mass of the top quark from pp collisions at 1.96 TeV observed 
with the Collider Detector at Fermilab (CDF) at the Fermilab Tevatron Run II. The events 
have the decay signature of pp — > ti in the lepton plus jets channel in which at least one jet is 
identified as coming from a secondary vertex and therefore a fe-hadron. The largest systematic 
uncertainty, the jet energy scale (JES), is convoluted with the statistical error using an in- 
situ measurement of the hadronic W boson mass. We calculate a likelihood for each event 
using leading-order ti and VK+jets cross-sections and parameterized parton showering. The 
final measured top quark mass and JES systematic is extracted from a joint likelihood of the 
product of individual event likelihoods. From 118 events observed in 680 pb~ of data, we 
measure a top quark mass of 174.09 ± 2.54 (stat+JES) ±1.35 (syst) GeV/c 2 . 



1 Introduction 

The top quark mass is a fundamental parameter in the Standard Model. Along with the mass of 
the W boson, the top quark mass provides the best indication for the value of the Higgs boson 
mass. This is a preliminary measurement of the top quark mass in pp collisions at y/s = 1.96 
TeV with the CDF detector^ at Run II of the Fermilab Tevatron. 

At the Tevatron, top quarks are mainly produced in pairs: 85% by qq annihilation and the 
rest by gg fusion. Top quarks decay before hadronization most of the time into a b quark and 
a W boson. The W bosons decay leptonically or hadronically and are used to topologically 
classify the final state of the ti system. We use events where one W boson decays hadronically 
and the other leptonically - the "lepton + jets" channel. Since taus have poor resolution, we 
include only electrons and muons and thus sample from 30% of the total branching ratio. 

This analysis uses a matrix-element analysis technique to efficiently extract information 
from the limited number of ti events. Each event enters the analysis with a weight derived 



from the differential cross-section for tt decay. We take into account all the possible jet-parton 
assignments and integrate over unknown quantities such as the momentum of the neutrino. The 
mass of the hadronically decaying W boson is measured and constrained with the current world 
average of 80.4 GeV/c 2 to measure the jet energy scale (JES) of the events. The jet energy 
scale is the largest source of systematic uncertainty in this analysis. The final top quark mass 
is extracted from a simultaneous fit of the top quark mass and JES. 

2 Data Sample &: Event Selection 

Events from the lepton + jets decay channel are selected requiring a single, high-transverse 
energy, well-isolated lepton; large missing transverse energy; and exactly four, central, high- 
transverse energy jets. (Two jets originate from the b quarks and two from the hadronically 
decaying W boson.) Of these jets, we require at least one to be "6-tagged," identified as orig- 
inating from a secondary vertex and thus the decay of a long lived b hadron. The primary 
vertex is the one from which tracks associated with the lepton emerge. The secondary vertex 
tag identifies tracks originating from a vertex displaced from the primary vertex and associated 
with a jet. To reduce the amount of non-VF background, we further require the leading jet and 
missing transverse energy not be collinear in the transverse plan for the lowest values of missing 
transverse energy passing our selection. Table H outlines this event selection. 



Table 1: Event Selection. 


lepton 


E T > 20 GeV (e), p T > 20 GeV/c (//) 


jets 


E T > 15 GeV, \r]\ < 2.0 


missing Et 


missing Et > 20 GeV 


6-tag 


> 1 jet coming from secondary vertex 


QCD veto 


0.5 < Acfx 2.5 (missing E T < 30 GeV) 



3 Method 

We write a likelihood for each event by combining a signal probability with a background 
probability. This method was first used to measure the top mass by the D0 collaboration 
during Run 1^. The likelihood is minimized for three variables: the top quark mass, the jet 
energy scale (JES), and the fraction of events consistent with our signal hypothesis (C s ), where 
x are the measured quantities: 

N 

C(Mto P , JES, C s ;x) oc Y[[C s P a (x; M top , JES) + (1 - C s )P w+iets (x ; JES)]. (1) 

i=l 

We first minimize the likelihood for C s with MINUIT and then perform a two dimensional fit to 
extract M top and JES. The signal probability, Pg, indicates how well an event describes leading 
order tt pair production and decay, and the background probability, -fV+jets> indicates how well 
an event describes the largest contributing background process, a leptonically decaying W boson 
plus extra jets. These probabilities are calculated by integrating over the parton differential 
cross-section, where the measured quantities are input to detector resolution functions, W(x, y), 
used to transfer the parton quantities, y, to measured quantities, and we also include parton 
distribution functions, /(%), for pp collisions: 



P(x) = — j da(y)W(x,y)f(q 1 )f(q 2 )dq 1 dq 2 . 



(2) 



We use an analytical form of the leading order matrix-element for qq — ► ti in the signal proba- 
bility and the sum of the W + 4 jets matrix-elements of the VECBOS Monte Carlo generator in 
the background probability. 

Our input measured quantities are the momentum of the lepton and the angles and energies 
of the jets. Regarding detector resolution, we consider the momentum of the electron or muon 
and the jet angles to be well measured. Thus, the detector resol ution in Equation [2] corresponds 
to the jet energy resolution, which is modeled from Monte Carlo^^ using a "transfer function," 
a mapping between jet energies and parton energies, Wj e t(x,y). Signal and background prob- 
abilities are summed over all possible permutations of jet and parton combinations, and the 
signal probability also considers different possible values of the transverse momentum of the ti 
system. 

We are sensitive to the jet energy scale through the mass of hadronically decaying W boson 
by constraining it with the world average W boson mass of 80.4 GeV/c 2 . We define the JES as 
a multiplicative scale factor applied to the energies of the two jets selected as the daughters of 
the W boson decay: 

Ejet = Ejef /JES. (3) 

The constraint comes from integrating over the event mass set by the daughters of the W boson 
and the JES in a Breit-Wigner using the world average value as the pole mass. We assume the 
JES determined for V7-jets also applies to b jets and assign a systematic uncertainty for the 
difference between the W and b jet energy scale. 

4 Systematic Uncertainties 



Table 2: Sources and Values of Systematic Error. 



Source of systematic uncertainty 


Magnitude (GeV/c 2 ) 


Residual jet energy scale 


0.42 


b-jet energy scale 


0.60 


Generator 


0.19 


Initial state radiation 


0.72 


Final state radiation 


0.76 


b-tag Et dependence 


0.31 


Background composition 


0.21 


Parton distribution functions 


0.12 


Monte Carlo statistics 


0.04 


Total 


1.35 



Table 121 lists the systematic uncertainties estimated from various Monte Carlo samples and 
re-weighting techniques. To first order, we fit out the JES systematic of our likelihood, but we 
also apply a residual systematic to cover an higher order effects, such as variations in the expected 
r\ or pt distribution. This higher order effect is estimated by shifting the input jet energies up 
and down one sigma as defined by the CDF Jet Energy and Resolution group^. The generator 
systematic uncertainty takes into account differences in the fragmentation and showering by 
comparing two different Monte Carlo models (PYTHIA and HERWIG). Possible biases originating 
from differences in the amount of initial- and final-state radiation between data and Monte Carlo 
are estimated using PYTHIA samples generator with lesser or greater amounts of radiation. The 
uncertainty on the parton distribution functions is evaluated as the sum in quadrature of the 
difference between MRST and CTEQ parton distribution functions; MRST with A QC d = 228 MeV 



CDF Run II preliminary (680 pb' 1 ) 

CO 
LU 

^> 1.05 
1 

0.95 
0.9 



165 170 175 180 185 

M top (GeV/c 2 ) 

Figure 1: Likelihood after minimization of C s parameter. 

and A.qcd = 300 MeV; and the variation of the 20 CTEQ6 eigenvectors. Systematic effects from 
the dependence of the 6-tagging on pr are evaluated by changing the dependence by one sigma. 
The background composition and modeling systematic error is the sum in quadrature of the 
largest variation when fluctuating the contribution of each individual background sample by 
100%; the change in the signal fraction by ±10%; and the largest variation due to changing the 
Q 2 scale in W + jet production. To understand the effects of limited statistics in the background 
sample, we divide the smallest sample (non-M^ background) in half and compare the results using 
each half separately. We repeat this procedure several times and histogram the difference. We 
take one-half the RMS of the distribution as a systematic error. 

5 Results 

The output likelihood is a simultaneous fit to three parameters: top quark mass, jet energy 
scale (JES) and signal fraction (C s ). We do not use any prior knowledge of JES or C s in the 
extraction of the top mass. The data used in the analysis corresponds to that collected in the 
period between March 2002 and September 2005, a total integrated luminosity of 680 pb -1 , 
where 118 events pass event selection. Figure ^ shows the fit to this data after minimization for 
C s as a function of Mt op and JES with different Aln£ contours. Our measurement is: 

M top = 174.09 ± 2.54 (stat + JES) ± 1.35 (syst) GeV/c 2 (4) 
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